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Single-molecule identification via electric current noise 
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Label-free and real-time single-molecule detection may aid the development of high-throughput 
biosensing platforms. Molecular fluctuations are a source of noise that often hinders single- 
molecule identification by obscuring the fine details of molecular identity. In this study, we 
report molecular identification through direct observation of quantum-fluctuation-induced 
inelastic noise in single organic molecules. We investigated current fluctuations flowingthrough 
a single molecule that is chemically connected to two electrodes. We found increased current 
oscillations synchronous to electric field excitations of characteristic molecular vibrational 
modes that contribute to inelastic electron tunnelling. This finding demonstrates a large 
contribution of charge interaction with nuclear dynamics on noise properties of single-molecule 
bridges and suggests a potential use of inelastic noise as a valuable molecular signature for 
single-molecule identification. 
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Recent advances in single-molecule detection techniques have 
revolutionized our understanding of biological phenomena 
and opened new opportunities for human gene analysis^ 
Unlike ensemble average measurements of macroscopic properties, 
molecular fluctuations generate a noise that often becomes a serious 
impediment for identifying important single-molecule signatures. 
In stark contrast, human sensory apparatus deals with the stochastic 
nature of biological systems by leveraging the unique ability to use 
noise for detecting weak electrical signals (via the stochastic reso- 
nance mechanism)^'^. Likewise, noise analyses have provided valu- 
able insights into single-molecule science. In particular, shot noise 
measurements have proven to be a powerful method for evaluating 
the Fano factors, as well as the number and transparency of chan- 
nels for charge transmission in individual molecules^'^^. In spite of 
these sophisticated experiments, little is known about noise charac- 
teristics in the inelastic regime, in which conducting electrons inter- 
act with molecular vibrational degrees of freedom. 

Here, we report the observation of inelastic contributions to cur- 
rent noise in single organic molecules connected to two metal elec- 
trodes. We analysed conductance fluctuations in a current-carrying 
molecular bridge under a high field. We found stepwise increases in 
the current oscillation amplitudes at distinct voltages. Our results 
suggest that this field-induced current noise originates from inelas- 
tic excitations of characteristic vibration modes inherent to target 
molecules, and thus manifests itself as a valuable molecular signa- 
ture for single-molecule identification. 

Results 

Current versus bias voltage in single-molecule junctions. We 

explored high-field conductance fluctuations in 1,6-hexanedithiol 
(HDT) single molecules trapped between two Au nanoelectrodes 
(HDT single-molecule junctions) using a self-breaking technique 
at 4K under a low applied voltage of 0.2 V (Fig. la,b; Supplementary 
Figs SI -S3). Figure Ic shows an I-V^ curve measured by a direct 
current (DC) method. The linear characteristics represent tunnelling 
electron transport through an HDT molecule, with a large highest 
occupied molecular orbital -lowest unoccupied molecular orbital 
(HOMO-LUMO) gap exceeding 6eV that yields a high carrier 
injection barrier at the metal-molecule interfaces when bonded 
to Au electrodes. In addition, we observe that current tends to 
become unstable with increasing V^. Correspondingly, the numerical 
differential conductance exhibits considerable fluctuations at high- 
regimes (Fig. Id). We note that the magnitude of the fluctuations is 
enhanced with increasing V^- 

Electric field dependence of current noise in HDT. To elucidate 
the underlying mechanism responsible for the electric-field-induced 
conductance fluctuations in single HDT molecules, we carried out 
consecutive 50 current measurements per condition at a sampling 
rate of 30 Hz during a bias sweep in a range of - 0.38 V<Vi,^ 0.38 V 
(Fig. le) and examined their average current <I> and current noise a 
(Fig. If). < / > is plotted as a function of in Figure 2a. We obtained 
linear <I>-V^ characteristics analogous to the /- curve shown in 
Figure Ic. Subtracting the elastic tunnelling current ^ from <I> yields 
the inelastic component A<I>= <J>-4^^'^1 We estimated ^ froni a 
linear fit to the </> - curve in a low bias range of - 0.005 V<Vb^ 
0.005 V; nonlinearity of elastic tunnelling contributions is anticipated 
to be negligibly small because of the relatively high tunnelling barrier 
at the electrode-molecule links expected for the insulating molecular 
wire of HDTs (Supplementary Fig. S4). By extrapolating linear I^-V^ 
to high biases (dotted line in Fig. 2a), we extracted A</> and plotted 
with respect to (Fig- 2b). The fact that |A< J>| rises above OA with 
increasing implies opening of active inelastic channels, rendering 
an increase in transmission of HDT molecules^ ^"^^ 

Meanwhile, a is found to increase with V^, manifesting more pro- 
nounced current fluctuations in the HDT single molecule under high 
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Figure 1 1 A measurement scheme for characterizing current fluctuations 
in single-molecule junctions, (a) A schematic illustration showing a 
measurement circuit used to investigate the high-field conductance 
fluctuations in Au-HDT-Au single-molecule junctions at 4.2 K. (b) lET 
spectrum of the HDT single-molecule junction formed by a self-breaking 
technique at 4.2 K. The peaks indicated by arrows can all be assigned 
to the lETS-active molecular vibrational modes (Supplementary Figs 
S1-S3 for details), (c) An example of l-V^ characteristics acquired 
for HDT single-molecule junctions at 1.3 mGo conductance state, (d) 
Differential conductance numerically derived from the l-V^ curve in b, 
revealing considerable fluctuations at high regimes, (e) 50-point 
current measurements conducted for evaluating the bias voltage-induced 
conductance fluctuations observed in HDT single-molecule junctions, 
(f) A magnified view of the plots in e. Black and red dots are the 50-point 
current data measured at each and the corresponding average 
current < / > . Error bars denote s.d. 



field (Fig. 2c). Moreover, although somewhat obscure at a low-V^ 
region, G-V^ plots demonstrate a staircase-like voltage dependence. 
We first note an overall linear response of a with respect to V^. This 
feature is generally interpreted by Hooge's empirical rule established 
for bulk materials that states a linear response of the resistance fluc- 
tuation-induced 1// noise amplitude to current^^'^^ The vaHdity of 
this rule in atomic- and molecular-scale conductors has been exam- 
ined recently for various nanomaterials, such as carbon nanotubes, 
graphenes and atom-sized contacts^^"^^ We subtracted the l//-like 
noise Gyj by linear fitting to C-V^ plots in the positive and nega- 
tive regimes. The resulting (a-o^/j^-V^ plots reveal clear steps 
at distinct (Fig. 2d). Measurements on another 40 HDT single- 
molecule bridges yielded similar results. Numerical differentiation 
as shown in Figure 2d yields a spectrum with pronounced peaks 
that emerged at positions corresponding to the a steps (Fig. 2e). 
These peaks are located at Vp = 30, 66, 127, 172, 300 and 337mV 
in the positive side and at Vp=-31, -79, -127, -199, -298 
and - 369 mV in the negative side (Table 1). The spectrum is there- 
fore almost symmetric with bias polarity. 

We examined the reproducibility of the characteristic a steps 
by statistically analysing G-V^ characteristics for 41 HDT single- 
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Figure 2 | Charge transmission and current fluctuations in a single- 
molecule tunnelling junction, (a) <I>-V^ characteristics obtained for 
an HDT single-molecule junction at 4 K. The dotted lines are the elastic 
current assumed to increase linearly with V^. (b) A< / > plots derived 
from </>-\/b curve in a. (c) o versus plots acquired for an HDT 
single-molecule junction at 4.2 K. (d) o plots after subtracting 1/f-like 
contributions o^f. (e) The di(G-Gyf)/diV^-V^ curve obtained numerically 
from d. The characteristic peaks are labelled numerically. 



Table 1 1 Peaks in the a-spectrum for an Au-HDT-Au single- 
molecule junction and an Au-Au single-atom chain. 



Peak no. 



a-spectrum 



mV 



Wavenumber 
(cm-^) 



Vibration 
modes 



Au-HDT-Au 



1 


30 


242 


220* 


v(Au-S) 


2 


66 


533 


650* 


v(C-S) 


3 


127 


1,025 


1,050* 


v(C-C) 


4 


172 


1,389 


1,300* 


ywCCH^) 


5 


300 


2,422 


2,575t 


v(S-H) 


6 


337 


2,721 


2,865* 


vs(CH2) 



Au-Au 
1 
2 



22 



65 
177 



81 1 
145t 



Transverse 
Longitudinal 



^Assigned molecular vibration mode frequency taken from ref. 21 (electron energy loss 
spectroscopy). 

i'Ref. 22 (Raman spectroscopy). 
:j:Ref. 30 (point-contact spectroscopy). 



molecule bridges at 1.3 mGo. A histogram constructed with data 
extracted from the a spectrum of single HDT molecules contains 
several peaks at almost symmetric positions with respect to zero 
bias (Fig. 3). As we see in Table 2, the statistical distribution fairly 
reproduces the a spectrum shown in Figure 2e. 




-0.3 -0.2 -0.1 0 0.1 0.2 0.3 
^p(V) 

Figure 3 | A peak voltage histogram constructed with data extracted 
from the a spectra of 41 HDT single-molecule junctions formed at 4.2 K. 

Peaks are numbered as 1-6 and V-6' in the positive and the negative side, 
respectively. 



Table 2 | Peaks in the histogram of Au-HDT-Au single- 
molecule junctions. 



Peak no. 



mV 



Au-HDT-Au 



1 


30 


242 


r 


-30 


242 


2 


78 


630 


r 


-80 


646 


3 


130 


1,049 


3' 


-134 


1,082 


4 


174 


1,405 


4' 


-199 


1,405 


5 


310 


2,503 


5' 


-318 


2,567 


6 


350 


2,825 


6' 


-358 


2,890 



Noise characteristics of gold single-atom contacts. It is interest- 
ing to note that the a peak profile agrees with what appears in the 
inelastic electron tunnelling (lET) spectrum^^ (Fig. lb). Indeed, 
we can assign a peaks to electron energy loss spectroscopy (EELS) 
and Raman-active vibrational modes^^'^^, as listed in Table 1. This 
implies that the high-field current fluctuations in single HDT 
molecules are of inelastic origin^^"^^ To verify the possible roles of 
inelastic excitations of vibration modes of single-molecule junc- 
tions on current fluctuations, we conducted a measurements on 
Au single-atom chains. Forming 1 Gq states using the self-breaking 
method, with a stretching speed of 6pms"^ at 4K, we exhibited 50- 
point current measurements in a bias range of -0.04 V<yb^ 0.04 V. 
The relatively small bias window was used because Au single-atom 
contacts tend to become increasingly unstable at bias voltage above 
0.05 V presumably because of the local heating and current-induced 
forces^l </>-Vi, plots in Figure 4a show linear characteristics 
reflecting the quasi-ballistic nature of electron transport through 
an Au single-atom contact^^. We extract the inelastic component by 
examining linear fitting to Figure 4a in a bias range from - 0.001 V 
to 0.001 V. In contrast to HDT molecules demonstrating increased 
transmission via active inelastic channels, A<I>-V^ plots in Figure 
4b reveal negative contributions of electron-phonon interactions in 
an Au single -atom chain on the electrical conductance^^ This dis- 
crepancy can be explained by considering the relative importance 
of the counteracting effects of elastic and inelastic processes in elec- 
tron transmission through single-atom and molecule junctions; 
whereas the inelastic mechanism leads to enhanced transmission. 
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Figure 4 | Inelastic contributions to electrical conductance and current 
fluctuations in a single-atom quasi-ballistic chain at 4.2 K. (a) <I>-V^ 
characteristics of an Au single-atom contact at 4 K. The dotted lines are 
the elastic component, (b) A< / > plots derived from the </> - V^ curve in 
a. (c) G-V^ plots measured for an Au single-atom contact formed at 4.2 K. 
(d) Plots of (cj-a^/f) with respect to and (e) the corresponding noise 
spectrum. Arrows point to the characteristic peaks in the spectrum at the 
positive range measured. 



elastic contribution involving virtual phonon excitations decreases 
transmission^^"^^. Theoretical studies have predicted a shift in the 
balance towards predominance of inelastic effects as transmission 
T becomes <0.5 and superior elastic effects for T becomes >0.5, 
the overall tendency of which agrees with the positive and negative 
A<I> acquired for HDT single-molecule junctions with T~10"^ and 
Au single-atom chains with T~l, respectively^^"^^. As for the high- 
field current fluctuations, a shows a stepwise manner analogous to 
that found for single HDT molecules (Fig. 4c,d). The corresponding 
a spectrum demonstrates peaks at VJ, = 8 and 22 mV (Fig. 4e). Sig- 
nificantly, we again find good accordance between these peak volt- 
ages and the characteristic vibration modes of single-atom wires at 
energies 10 and 18meV reported in previous scanning tunnelling 
microscopy studies^^'^^ (Table 1). It is therefore corroborated that 
bias voltage-induced conductance fluctuations can serve as a unique 
molecular signature for single-molecule identifications. 

Discussion 

The above results strongly indicate that bias voltage-dependent cur- 
rent oscillations are phonon- assisted phenomena relevant to electron 
vibration interactions in atomic/molecular bridges. Inelastic tunnel- 
ling in molecular systems is known to cause a local heating^^'^^"^^. 
A certain amount of heat is anticipated to pile up in the molecular 
junctions every time a new inelastic channel sets in. Although this 
excess heat dissipates via lattice conduction through two electrodes 
connecting a molecule, local heating effectively raises the junction 
temperature T^^^ above the bath temperature^^'^^"^^. This may yield a 
stepwise increase in the thermal noise with increasing V^, which 



explains a staircase-like increase of a observed in single HDT mol- 
ecules and also in Au single-atom wires. It has been theoretically pro- 
posed that the heating rate can be deduced by analysing the bias voltage 
dependence of Johnson-Nyquist noise in molecular junctions^^'^^. We 
therefore sought to extract information regarding electrical heating 
in HDT junctions from the a plots in Figure 2d. We evaluated T^^f 
by tentatively considering Gt = {c-GyX where is the thermal noise 
component of a, and assuming a^^^ v T^g. The relative increase in Oj at 
Vb = 0.4 V from the zero-bias value is about 20, from which we obtain 
Tgff- 1,600 K at a bath temperature of 4K. The estimated T^g is obvi- 
ously too high and contradicts the fact that an HDT molecule can be 
held stable during a bias sweep without being burned out. 

The overestimation of T^^^, which is also marked by the fact that the 
lET spectrum can be observed, suggests significant inelastic noise 
in the G-G^f of single HDT molecules derived from the stochas- 
tic nature of phonon-assisted electron tunnelling and concomitant 
fluctuations in the phonon-induced energy levels^^'^^'^^. Shot noise 
measurements have proven to be a powerful method for evaluat- 
ing Fano factors, as well as the number and transparency of charge 
transmission channels in atomic/molecular conductors^'^*^'^^. A theo- 
retical model predicts a linear increase in the power spectral density 
of this shot noise with the bias voltage in a quasi-ballistic Au single- 
atom chain under a condition T^^> hv/k^, where V is the frequency 
of phonon modes^^. Indeed, we found a stepwise linear increase in 
G-Gyf after each vibration mode excitation above 0.01 V, as shown 
in Figure 2e, which in turn suggests a significant local heating in an 
Au chain that raises its T^^ to > K at = 0-02 V by virtue 

of the poor contact-cooling efficiency via lattice heat conduction at 
4 K. Furthermore, the positive contributions of inelastic channels to 
noise amplitude found in the present study suggest a weak electron- 
phonon interaction and symmetric molecule-metal electronic cou- 
pling in HDT single-molecule junctions^^. (It should be noted that 
our discussions on inelastic noise is based on a fixed sampling rate. 
Although the noise amplitude is anticipated to change at different 
data acquisition speeds, the qualitative arguments on the voltage- 
specific change in noise amplitudes attributable to the electron- 
vibration interactions in atomic/molecular conductors still hold.) 
These results indicate that vibrational mode-specific inelastic noise 
in single molecules can be more than a useful molecular signature 
to identify single-molecule junctions: with elaborate theoretical 
models, it can serve as a unique probe for detailed investigation of 
the electronic properties of single-molecule bridges. 

Methods 

Fabrication of mechanically controllable break junctions (MCBJs). Lithographi- 
cally defined nano-MCBJs were exploited to form single-molecule bridges. Fabrica- 
tion procedures of nano-MCBJs are described elsewhere^^. In short, we spin-coated a 
polyimide on a phosphor bronze substrate for electrical insulation. Sub -micrometer- 
scale Au junctions were then patterned using an electron beam lithography tech- 
nique and subsequent lift-off processes. The sample was then exposed to reactive 
ion etching to remove polyimide underneath the junctions. As a result, we obtained 
Au nanobridges with free-standing length of about 2 |im. In experiments, we first 
bent a substrate with a three-point bending mechanism and ruptured Au junctions 
immersed in a dilute toluene solution of test molecules. The sample geometries were 
designed to provide the damping factor r=Ad/AL-3xlO~'^, where Ad and AL denote 
the displacements of the MCBJ electrodes and the substrate deflection, respec- 
tively^*^. Here, we used HDTs, benchmark molecules the electron transport proper- 
ties of which have been studied intensively^"^'^". Subsequently, we evacuated the 
sample chamber to remove the solvent, leaving HDT molecules adsorbed on a fresh 
fracture surface of Au junctions. When the vacuum level reached below 10"^ torr, we 
performed opening and closing of Au junctions repeatedly for hundreds of cycles at 
room temperatures. At the same time, we examined the evolution of Au-HDT-Au 
structures by monitoring temporal conductance changes of junctions under a junc- 
tion-stretching speed of 0.6 pms"^ using a picoammeter (Keithley 6487) under 
a DC bias voltage of 0.2 V. The measured G-t curves revealed conductance steps 
at < IGq, suggestive of formations of Au-HDT-Au systems (Supplementary Fig. SI). 
After this, we slowly cooled the junctions to 4K. 

Formations of Au-HDT-Au single-molecule bridges. We created an HDT 
single- molecule bridge using a self-breaking method^ under Vj = 6pms"^ and 
V^ = 0.2y at 4K. Conductance decreased in a stepwise manner. We stopped junc- 
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tion stretching when a well-defined G-plateau at 1.3x10"^ Gq was detected right 
after a conductance drop from 1 Gq, which indicates opening of an atom- sized 
electrode gap (Supplementary Fig. S2). 

Inelastic electron tunnelling spectroscopy measurements. To verify the 
presence of HDT molecules between two electrodes, we exhibited inelastic electron 
tunnelling spectroscopy (lETS) analyses by measuring dl/dV^ in a bias range of 
OAy<V^<OAy using a lock- in amplifier (SR830, Stanford Research Systems) with 
a reference AC voltage of 10 mV^^^ at 1 kHz after forming and holding the single- 
molecule conductance state (Supplementary Fig. S3). The differential conductance 
was found to increase in a stepwise manner with V^, indicating the opening of 
inelastic channels. d^I/dV^^-V^ characteristics were obtained numerically from the 
measured dl/dV^ curves. Pronounced peaks were observed in the lET spectrum 
that can be assigned to the lETS-active vibrational modes of Au-alkanedithiol-Au 
bridges. 
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